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Abstract 

We give a new derivation of the quasinormal frequencies of Schwarz- 
schild black holes in d > 4 and Reissner-Nordstr0m black holes in d = 4, 
in the limit of infinite damping. For Schwarzschild in d > 4 we find that 
the asymptotic real part is ^Hawking l°g(3) for scalar perturbations and 
for some gravitational perturbations; this confirms a result previously 
obtained by other means in the case d = 4. For Reissner-Nordstr0m in 
d = 4 we find a specific generally aperiodic behavior for the quasinor- 
mal frequencies, both for scalar perturbations and for electromagnetic- 
gravitational perturbations. The formulae are obtained by studying the 
monodromy of the perturbation analytically continued to the complex 
plane; the analysis depends essentially on the behavior of the potential 
in the "unphysical" region near the black hole singularity. 
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1 Introduction: why do we study 
quasinormal modes? 

Recently it has been proposed [1] that the asymptotic behavior of the high 
overtone black hole quasinormal frequencies captures important information 
about the spectrum of black hole observables and about quantum gravity 
in general. More precisely, it has been suggested that the asymptotic real 
part of the black hole quasinormal frequencies coincides with the frequency 
emitted by a black hole whose area decreases by an amount that is natu- 
ral from the point of view of discrete approaches to quantization of gravity 
such as loop quantum gravity (LQG). While such approaches remain highly 
speculative, it is nevertheless interesting to understand how far this argu- 
ment can be pushed; moreover, even if this particular line of inquiry has to 
be abandoned, the asymptotic spectrum of black hole excitations is a fun- 
damental property of the black hole which should eventually be good for 
something. This paper is an attempt to simplify and extend the calculations 
of the relevant "experimental data." 

Let us briefly review the recent developments in loop quantum gravity 
which have led us to investigate this asymptotic spectrum. (See [2, 3, 4] for 
increasingly lengthy reviews of LQG and [5, 6] for the original derivations of 
the area quantization law in LQG.) In LQG, the natural quantum Aq of the 
horizon's area is proportional to the so-called Barbero-Immirzi parameter 7 
[7, 8], a parameter in LQG which is arbitrary from the point of view of the 
microscopic definition of the theory, but whose value must be chosen in such 
a way that the black hole entropy comes out correctly. This comes about as 
follows: one imagines that each quantum of area Aq carries a piece of discrete 
information. It can be found in k different microstates, and therefore it 
carries entropy log A;. This must be equal to the known contribution of this 
area quantum to the black hole entropy, namely 1 Aq/4G, so we find that 

A = AGlogk. (1) 

LQG predicts that k = 2j( min ) + 1, where j( mm ) is the minimal allowed 
spin carried by the spin networks. This means that At = 2 if the basic 
gauge group of LQG is taken to be SU(2), or k = 3 if the gauge group is 
SO(3). However, from the LQG point of view there seems to be no reason 
to choose this strange number log(2) or log(3) (more precisely, log(2)/7r\/3 
or log(3)/27T\/2) for the Barbero-Immirzi parameter. 

On the other hand, seemingly unrelated, one can consider the "quasinor- 
mal modes" of the black hole background. Quasinormal modes are pertur- 
1 In this paper, we use particle physics units with c = h — 1. 
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bations of the background which are allowed to have complex frequency but 
whose boundary conditions must be "purely outgoing" both at the horizon 
and at infinity; we will review their precise definition in Section 2. This 
boundary condition singles out discrete complex frequencies u n . As n — > oo, 
the quasinormal frequencies of gravitational perturbations are known to be- 
have as 

w« = T Ha wking(27rz(n - 1/2) + log(3)). (2) 

(Varying from the usual convention in general relativity, we parameterize 
the frequencies in terms of ^Hawking) which is equal to 1/SttMG for a four- 
dimensional Schwarzschild black hole.) The imaginary part of (2) is not 
hard to understand. The quasinormal modes determine the position of poles 
of a Green's function on the given background, and the Euclidean black hole 
solution converges to a thermal circle at infinity with the inverse temperature 
(5 = 1/Tnawking; therefore it is not too surprising that the spacing of the poles 
in (2) coincides with the spacing 27riTH aw king expected for a thermal Green's 
function. 

However, we are not aware of a similar classical explanation of the term 
log(3) in (2). This term is remarkable because, as was first observed in [1], 
it is precisely the same strange number one needs to occur in the Barbero- 
Immirzi parameter of LQG if the gauge group is taken to be SO (3). There 
are also some heuristic arguments which suggest that there might indeed be 
a reason to relate the two, identifying the quasinormal mode with a process 
in which the black hole ejects a single spin-network edge and thus reduces 
its area by A$. 

The result (2) was originally obtained by numerical calculations ten years 
ago [9, 10, 11, 12]. At that time the term log(3) was only known to an 
accuracy of six decimal places. In 1998, Hod conjectured [13] that this 
number was exactly log(3). His conjecture was recently proven in [14], using 
the methods of continued fractions initiated by Leaver [9] and further refined 
by Nollert [10]. 

In this paper we would like to reproduce the results of [14], using sim- 
pler functional methods which also allow us to extend the analysis to the 
cases of charged and higher-dimensional black holes. We will find that the 
answer ^Hawking log(3) is in fact universal for Schwarzschild black holes in 
all dimensions, at least for scalar and some gravitational perturbations. For 
Reissner-Nordstr0m in d = 4 we will find a more complicated generically 
aperiodic behavior, both for scalar and electromagnetic-gravitational per- 
turbations. 

Our computation depends on an analytic continuation in the radial coor- 
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dinate r, so it is worth mentioning that analytic continuation in r has been 
applied to the determination of quasinormal modes many times before; in 
particular, it was the basis of the numerical calculations in [11] which con- 
firmed that Reoj — > Tnawking log(3) for Schwarzschild. A review of the use of 
complex coordinate methods to compute black hole quasinormal frequencies 
has recently appeared in [15], to which we refer the reader for more extensive 
references. We note that in our calculation of the leading-order asymptotic 
frequencies no phase integral technology will be required; it will be enough 
to use simply the first-order WKB approximation. 

We believe that this puzzling behavior of the Green's functions at large 
imaginary frequencies should contain some important information about 
black holes and perhaps even quantum gravity. 

2 Quasinormal frequencies for d = 4 Schwarzschild 
black holes 

2.1 The differential equation 

In this section we will be concerned with properties of certain perturbations 
of a four-dimensional Schwarzschild geometry, described by solutions of the 
differential equation [16] 

+ V(r) - uo^j 4>{r) = (3) 
on the interval 1 < r < oo, with the "Regge- Wheeler" potential 

v<0-(.-i)(^ + i^). 

The equation (3) describes the radial dependence of a spin-j perturbation 
of the background, with orbital angular momentum I and time dependence 
e~ lu>t . We have set 2GM = 1; in these units the horizon is at r = 1. 

The numerator (1 — j 2 ) which appears in (4) characterizes the spin j of 
the perturbation for j = 0, 1,2, i.e. for perturbations by scalar or electro- 
magnetic test fields as well as perturbations of the metric itself. (Actually, 
in the cases j = 1,2, one has to expand in tensor spherical harmonics and 
then (4) is only correct if we consider an "axial" perturbation; see e.g. [18].) 
It is nevertheless interesting to study the problem for general j, because 
our calculation relates this coefficient to Bessel's equation and shows that it 
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is very natural to write the numerator in this form. It might be useful to 
summarize its values for the most important examples of j: 



scalar perturbation j ' = 
electromagnetic perturbation j = 1 (5) 
gravitational perturbation j = 2 



The factor (l — ^) which occurs in (3), (4) is the "warp factor" in the 
Schwarzschild geometry, which in Schwarzschild coordinates is simply the 
time component —goo of the metric. 

We wish to study (3) on the physical region 1 < r < oo and we begin 
by describing the asymptotic solutions. To simplify the differential operator 
which appears it is convenient to define the "tortoise coordinate" x with the 
property that 

l^- 1 



dx = - -j dr. (6) 
Integrating (6) gives 

x = r + log(r- 1). (7) 
Then x takes values — oo < x < oo, and (3) becomes 



2 



^2+V[r{x)]-u 2 )i- = 0. (8) 



Since V[r(x)] — ► as x — > ±oo, the solutions behave as tp ~ e ±luJX at infinity. 



2.2 Analytic continuation 

In our analysis we will find that it is essential to extend (3) beyond the usual 
physical region 1 < r < oo. Continued to the whole complex r-plane, (3) 
is an ordinary differential equation with regular singular points at r = 0, 
r = 1 and an irregular singular point at r = oo; and by the general theory of 
differential equations, any solution of (3) in the physical region extends to a 
solution on the r-plane. However, this solution may be multivalued around 
the singular points r = 0, r = 1. This multivaluedness will play a crucial 
role in our analysis, which is essentially concerned with the computation 
of the monodromy around a particular closed contour in the r-plane. For 
convenience, to avoid having to deal with multivalued functions, we can put 
branch cuts in the r-plane from r = and r = 1 and require that ip(r) 
satisfy (3) everywhere except at the cuts. In this approach the monodromy 
is defined by the discontinuity across the cut. 
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Figure 1: The r-plane, with regions 
where Re x < denoted by a darker 
color, and a convenient choice of 
branch cuts shown. 



v. 



We will find it very convenient to exploit the simpler representation (8) 
for the differential equation. After passing to complex r this introduces an 
extra complication, since the relation (7) between x and r is itself multivalued 
because of the term log(r — 1). As a result, as we travel around a contour in 
the r-plane, the corresponding contour in the x-plane may not be closed. In 
the language of branch cuts we would say that x jumps when r crosses the 
branch cut which emanates from r = 1. 

Although x is not uniquely defined as a function of r, the ambiguity is 
only in the imaginary part, so Rex still makes sense. It will be important 
for us to know the sign of Rex in the various regions of the r-plane; the 
structure is as shown in Figure 1. We also show a convenient choice of branch 
cuts (to economize, one should take the branch cut to be the same for ip(r) 
and for log(r — 1)) although nothing we say will depend on this choice. 

2.3 Definition of quasinormal modes 

Because the potential V(r) is positive and decays to zero at both ends, there 
are no discrete normalizable bound states. Nevertheless we may look for 
discrete quasinormal states, counterparts of quasistationary states in quan- 
tum mechanics because their frequency is also allowed to be complex. The 
mathematical theory of these quasinormal modes is more intricate than for 
ordinary normal modes; see e.g. [20] for a more precise discussion than we 
will give here, as well as an excellent review of the substantial literature on 
quasinormal modes of black holes. 

Roughly, a quasinormal mode with frequency u is supposed to be a solu- 
tion of (8) where we require "outgoing" behavior at both ends of the physical 
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region: 




as 



as x — > +00, 



x — > —00. 



(9) 
(10) 



Since the differential operator in (3) would be self-adjoint acting on normal- 
izable wavefunctions, and hence would have to have real eigenvalues, there 
are no quasinormal modes with Imu < 0. So we take Iuilo > 0, and then 
in (9), (10) we are asking for the exponentially decreasing component to be 
absent at both ends. It is actually quite delicate to define exactly what this 
means! If V[r(x)] were strictly vanishing at x — > ±00 then one could just 
require that tp be strictly proportional to e ±tuJX as appropriate, but for more 
general V[r(a;)] we do not know any way of defining the boundary condition 
using only the behavior of ip on the real line. This problem was also dis- 
cussed in [19] where it was suggested that the appropriate condition would 
be obtained by finding an analytic extension of ip to Imw < and demand- 
ing exponential decrease there; this was further shown to be equivalent to 
looking for poles in an appropriate Green's function. In [9] another criterion 
was given in terms of the coefficients of a power series representation for ip; 
this criterion is what was used in the analytical calculation of the asymptotic 
quasinormal frequencies in [14]. 

Here we wish to use another way of defining the boundary condition, by 
analytic continuation in r instead of lo. (The method we describe below has 
been studied before, e.g. in [21]; see also [15] for a recent discussion and 
additional references.) By considering complex r we can easily define the 
boundary condition at the regular singular point r = 1, because the functions 
e ±iujx are distinguished by their monodromy there. Namely, e ±lwx = (r — 
l) ±la; e ±la;r , so to get tp ~ e tuJX near r = 1 we simply require that ip{r) have 
monodromy e _27rw on a counterclockwise circle around r = 1. What is more 
subtle is to define the boundary condition at r = 00, but this can be done 
by analytically continuing ip(r) via "Wick rotation" to the line Imcox = 0, 
where the asymptotic behavior of the solutions is purely oscillatory and hence 
one can really pick out a particular solution by specifying its asymptotics. 
In the limit where to is almost pure imaginary, the line Imwi = is slightly 
tilted off the line Rex = 0. There are two possible directions for the Wick 
rotation and we choose the rotation through an angle smaller than tt/2. For 
the time being we assume Re a; > 0; in this case x = 00 gets rotated to 
lux = 00, so our boundary condition "at r = 00" is actually 



V>(r) ~ e %wx as ujx — > 00. 



(11) 
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2.4 Computation of quasinormal frequencies 

Having defined the quasinormal modes we now proceed to study their fre- 
quencies, in the limit Imw > Rew. Indeed, until the end of this computa- 
tion we will behave as if uj were purely imaginary. 

We will compute in two different ways the monodromy of if)(r) around 
the contour shown in Figure 2, beginning and ending at the point marked 



Figure 2: Contour for the calculation of quasinormal frequencies, 
shown in terms of r and in terms of z — x — iir. The singularities 
at r = (cquivalently z = 0) and r — 1 (equivalently z = — oo) are 
marked. The rotation by 3tt/2 in the r plane becomes a rotation by 
37r in the z plane as follows from (12). The monodromy of log(r — 1) 
around r — 1 causes a 2m discontinuity in the path in the z-planc. 
The regions with Rex < are denoted by a darker color. The wave 
function is well approximated by a sum of Bessel functions in the blue 
circle around r = 0. 

It will prove convenient to use the coordinate z = x — m instead of 
x. Actually, this definition is slightly imprecise since we have not specified 
which branch of the logarithm in x(r) we should choose. Depending on the 
branch, the contour in the z-plane will approach some point z = 2mn as r 
approaches 0; for convenience we choose the branch so that n = 0. 

First we compute the monodromy by matching the asymptotics along 
the line Imwx = 0. Beginning at A, the plane wave behavior (11) can be 
extrapolated toward the interior since the term uj 2 dominates V in (8) away 
from a small neighborhood around r = 0. As we extrapolate up from A we 
reach r = at z = 0, and near this point we have 



A. 




z = r + log(r — 1) — iir ~ r , 



(12) 
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so in a neighborhood of r = we may write 

V[r(z)] ~ ( 13 ) 
Thus (3) may be simplified near r = to 

("l? + ^-" 2 )^> = ' < 14 > 

This equation can be exactly solved in terms of Bessel functions; for generic 
j, the general solution is 

ip(z) = A + c + \fuzJj/2(wz) + A^C-\/ r uJzJ_j/ 2 (^z). (15) 
We choose the normalization factors c + , c_ so that the asymptotics are 

c±y/uzJ±j/2(ujz) ~ 2cos(u>z — a±) as ujz — > oo, (16) 
where we defined 

a± = j(l±j). (17) 

The phase shifts by a± will be very important for our analysis. Using (15) 
and (16) together with the boundary condition (11) which states that e lwz 
must be absent as ^ oo, we have 

A +e -' ia + + A_e~ ia - = 0, (18) 

while for the coefficient of e~ %u)Z we have 

ijj(z) ~ (A+e ia + + A_e ia -) e~ iwz as ujz -► oo. (19) 

To follow the contour through to B we now have to turn through an angle 
37r/2 around r = 0, or equivalently 37r around z = 0. Using the fact that 

J± j/ 2(uz) = z ± f' 2 <l>(z), (20) 

where is an even holomorphic function of z, after the 3n rotation the 
asymptotics are 

c±y/ujzJ±j/ 2 {wz) ~ e 6ta± 2 cos(— ujz — a±) as ujz — > — oo. (21) 

Then using (15), (21) we get the asymptotics at B, 

i/>(z) ~ (A+e 5ia + + ^_e 5ia -) e - *"* + (A + e 7ia + + A_e 7ia -) e icJ2 (22) 

as ujz — > — oo. Finally, we can continue along the large semicircle in the 
right half-plane from B back to A. Since we are in the region where uj 2 
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dominates the potential in (8) we essentially just have plane wave behavior; 
this approximation can be trusted at least for the coefficient of e~ tuJZ , which 
therefore remains unchanged as we continue back to A and complete the 
contour (however, the same cannot be said for the coefficient of e tuJZ , which 
makes only an exponentially small contribution to V( z ) i n the- right half- 
plane and hence can be modified by the small correction terms in the WKB 
approximation.) In the end we find that the monodromy around the contour 
multiplies the coefficient of e~ tuJZ by a factor 

A + e 5ia + + A_e 5ia - e 6ia + - e 6ia - sin 37^/2 . . . . 

A + e ta + + A-e ia - e zux + — e zta - sm7r//2 



On the other hand, we can also obtain the monodromy action by observ- 
ing that the only singularity of ip(r) or e~ Vj}Z inside the contour is at r = 1. 
By our boundary condition going around r = 1 clockwise just multiplies 
ip(r) by e 27rw , but it multiplies e~ luJZ by e~ 27ra; , so the coefficient of e~ wz in 
the asymptotics of ip(r) must be multiplied by e AnuJ after the full round trip. 
Comparing with (23) gives 

e 47ru, = _(l + 2cOS7Tj). (24) 



Note that if we had chosen Re to < we would have had to Wick rotate 
in the opposite direction, hence set our boundary conditions at B instead of 
A, and run the contour in the opposite direction; hence we would have gotten 
(24) with e^ 4 ™ instead of e 4 ™. This is as expected since if u is a quasinormal 
frequency then —to* must also be one; the corresponding wave functions are 
simply complex conjugates of one another. The two possibilities can be 
summarized by the equation 

exp (e (Re uj) Attuj) = -(1 + 2 cos irj), (25) 

implying 

|4vrReo;| = log|l + 2cosvrj|. (26) 

In the case j = 0, 2 we recover the numerical result of [9, 10, 11, 12] which 
was so important for [1, 13]: 

4ttuj = (2n + l)vri±log(3). (27) 

If our analysis of the dependence on e(Rew) is correct, it implies |1 + 
2 cos irj\ > 1. Hence for j £ (1/2, 3/2) + 2Z, there cannot be any asymptotic 
quasinormal modes, with the possible exception of the case j 6 2Z + 1. (This 
case is exceptional in our analysis for another reason as well: if j = 1, then 
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our form (13) for the leading singularity vanishes. A more careful analysis 
would be desirable to check whether our results still hold in this case.) 

We should mention that (25) is equivalent to a more complicated equation 
derived in [14], namely 

tan(7ricj) tan(-7uu; + irj/2) tan(7rzo; — irj/2) = ie(Reuj) (28) 

This equivalence can be shown by multiplying (28) by 

cos(iriuj) cos(7riu; + 7rj/2) cos(7riu; — 717/2) (29) 

and expressing the sines and cosines in terms of exponentials. 

2.5 Some convenient notation 

For later convenience we now describe the preceding computation in a more 
compact notation. First write 

{m) A = (A + e ima+ + A.e ima - ) (30) 

and 

(m, n) A = (m) A e-^ 2 + (n) A e^ 2 . (31) 

Then the asymptotic matching through the Bessel region which we did by 
comparing (16) and (21) shows that a Sir clockwise rotation in z transforms 
the asymptotic behavior by 

3vr : (m,n) -> (n + 6,m + 6), (32) 

where the swap between m and n arises because of the swap ujz <-> —ojz 
needed to compare (16) and (21). 

In the case of four-dimensional Schwarzschild we read off the initial 
asymptotics from (16) as tp ~ (1,— 1)a- Since the coefficient of e wx had 
to vanish at A we got the constraint 

(-1)a = 0. (33) 

The 37r rotation (32) gave 

(1,-1)a-(5,7)a, (34) 

so following the coefficient of e~ tuJX back around the semicircle gives its mon- 
odromy as (5)a/(1)a, which we computed using the constraint (33), obtain- 
ing 

= -(l + 2cosvrj). (35) 
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3 Quasinormal frequencies for other black holes 

Having reproduced the result of [14] for d = 4 Schwarzschild, we now apply 
the same technology to two other classes of black hole. The generalizations 
are surprisingly straightforward. 



3.1 Higher-dimensional Schwarzschild black holes 

First we consider Schwarzschild black holes in arbitrary d > 4. We will verify 
a conjecture from [14]: namely, the asymptotic quasinormal frequencies of 
the Schwarzschild black hole with Re a; > satisfy 



(36) 



where (3 is the Hawking temperature of the black hole. (Strictly speaking, we 
will prove this only for scalar perturbations and for certain tensor modes; we 
do not know the appropriate potentials for general perturbations in d > 4.) 

To write the generalization of (3) to d > 4 we need the warp factor for 
d-dimensional Schwarzschild, 



-5oo = f(r) = 1 - ^3, 



(37) 



where we again chose units in which the horizon is at r = 1. Then the 
perturbation satisfies 



fir, 



d_ 

dr 



V(r) 



V>(r) = 0, 



where V(r) is known for scalar perturbations to be [22] 
V(r) = f(r) 



l(l + d-3) + (rf-2)(d-4)/(r) + (d - 2)/'(r) 



4r 2 2r 



As in the d = 4 case, this V(r) grows large only at r 
computation gives the coefficient of the singularity as 



(38) 



(39) 



0, and direct 



V(r) ~ - 



(d-2) 



4r 



2d-4 ' 



(40) 
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Figure 3: Regions with Re z < 
(darker) in the r-plane for a 
six-dimensional Schwarzschild black 
hole, together with a contour for the 
calculation of quasinormal frequen- 
cies. The left part of the r-plane is 
shown in shadow because Re z can 
only be uniquely defined in the right 
part; the extra fictitious horizons are 
marked by crosses. 



Now we want to describe the tortoise coordinate in d spacetime dimen- 
sions. Prom the warp factor (37) we get the requirement 



dz=\l 
which can be integrated to give 



1 



,-d— 3 



-1 



dr, 



(41) 



^ e 2nij/(d-3) 

r + Z-> J o 1 °g( 1 

j=0 



d-3 



re 



-2itij/(d-3) > 



(42) 



where we chose the additive constant so that z = for r = 0. From (42) 
we obtain the behavior of Rez as shown in Figure 3. Unlike the situation 
in d = 4, (42) includes phases multiplying the logarithms which come from 
the fictitious "horizons" r = e 2m i/( d ~ 3 \ so that Rez is not well defined for 
arbitrary r; but so long as we stay in the right part of the r-plane it is well 
defined. 



Near r = 0, d — 2 terms in the Taylor expansion of (42) near r = cancel 



and 



„d-2 



>~— r (43) 

generalizing the behavior (12) in the d = 4 case. Actually, this cancellation 
can be seen directly from (41), by neglecting the term 1 on the right hand 
side. 

Now we can express the singularity of the potential in terms of z; com- 
bining (40) and (43) we find simply 



V[r{z)\ 



1 



(44) 



Note that this is exactly the same as the behavior for scalar perturbations in 
d = 4; namely, it agrees with (13) at j = 0. The other case for which we know 
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the potential is a particular family of tensor metric perturbations studied in 
[23], which exist only in d > 4; for these the leading singularity corresponds 
again to j = 0. (One might have expected j = 2, but in our final result j 
appears only in cos7rj, so that j = and j = 2 are indistinguishable.) We 
conjecture that with some suitable interpretation this form for the potential 
will actually persist for all j, so that the —1 may be replaced by j 2 — 1 in 
(44), but the skeptical reader may just set j = in the rest of this section. 

At any rate, writing the potential in this form the Bessel function behav- 
ior of the solutions near r = 0, expressed in terms of z, is just as it was in 
d = 4, and again we have plane wave behavior away from a neighborhood of 
r = 0. Furthermore, although the left part of the r-plane is rather compli- 
cated, the behavior of Rex on the right part is essentially unchanged from 
d = 4. So we can repeat the argument from the previous section. The cal- 
culation of the monodromy by matching asymptotics in the z-plane works 
exactly as in four dimensions, while the factor e 4nuJ from the monodromy 
around r = 1 is modified to e 47rw /( d - 3 ) ; because of the coefficient l/(d — 3) 
multiplying log(l — r) in (42). 

Fortunately, this same factor appears in the formula for the Hawking 
temperature of a d-dimensional black hole. Namely, recall that the Hawking 
temperature can be computed as the surface gravity of the black hole: 

^Hawking = ^ = (45) 

giving (3 = Amj{d — 3). Hence the monodromy of the coefficient of e~ %wx 
around the contour is e 13 ^, and comparing this to the result from Bessel 
asymptotics gives 

e#" = -(1 + 2cos7T7), (46) 

as desired. 

In fact the agreement between the surface gravity and the coefficient 
of the logarithm in the tortoise coordinate is not a coincidence and holds 
more generally. This is because a horizon is characterized by the condition 
/(rff) = 0, and the tortoise coordinate near ru then obeys 

dz = / 7^ ~ / t r 77? — r = -jf, — r log(r - r H ) (47) 

J f(r) J (r-r H )f'(r H ) f'(r H ) 

while the surface gravity is / / (r//)/47r. 

Our result implies that for scalar and some gravitational perturbations 



Reuj = ±T Ha wkinglog(3). 



(48) 
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In the revised version of [24], considerations of Bohr-Sommerfeld quantiza- 
tion (as originally proposed by Bekenstein and amplified e.g. in [25]) led 
to the conjecture that Rew = ±7H aw king log(fc) for some k; our calculation 
establishes this conjecture and furthermore determines the integer k = 3. 



3.2 Four-dimensional non-extremal Reissner-Nordstr0m black 
holes 



Finally we want to extend our results to non-extremal Reissner-Nordstr0m 
black holes in d = 4. In this case the warp factor is 

/M-l-^ + £=<**> (49) 

where we have fixed our units by setting GM = (1 + k)/2, Q 2 = k, so that 
the outer horizon is at r = 1 and the inner horizon is at r = k < 1. The 
tortoise coordinate is therefore determined by 

dZ = dr (r-l)(r-kY (50) 

which integrates to 

. log(l - r) - /c 2 log(l - r/k) 

z = r + — (51) 

where we have again fixed the constant so that z = at r = 0. Then near 
r = we find 

2 ~5- < 52 > 

The form of the potential can be found in [29]: both for scalar perturbations 
and for axial electromagnetic-gravitational perturbations (which are coupled 
in this background, so we cannot talk about electromagnetic or gravitational 
perturbations separately) the potential is singular only near r = 0, and the 
leading singularity is 

, . J —2k 2 1 V 6 scalar perturbation . . 

| +4k 2 /r 6 electromagnetic-gravitational perturbation 

Using (52) we can rewrite this in the familiar form 

nr(z)] ~ ^ ( 54 ) 
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(with all the k dependence disappearing) and, remarkably, the values of j 
which come out are nice: 



Figure 4: Regions with Re (z) < (darker) in the r-plane for 
four-dimensional non-extremal Reissner-Nordstr0m, with a contour for 
the calculation of quasinormal frequencies, shown both in the r-plane 
and the z-plane. 

Hence the Bessel function behavior near r = will be unchanged from the 
previous examples, and our calculation will proceed in an almost identical 
manner. However, in this case we must use a slightly more complicated 
contour, going around a "lobe" as shown in Figure 4; the point is that we 
want to capture only the outer horizon at r = 1, where we know the boundary 
conditions and hence know the monodromy, and avoid the inner horizon at 
r = k. As a result this computation is only valid in the non-extremal case 



The corresponding picture in the z-plane looks rather complicated be- 
cause a given z may or may not lie in a "Bessel region" depending on whether 
or not the r-plane contour is passing near r = (the potential is a multi- 
valued function of z but an ordinary function of r.) In Figure 4 this is 
indicated by letting the contour pass over or under the blue circles when z 
is not passing through a Bessel region; it is easiest to understand that figure 
by redrawing it yourself, following the given r-plane contour starting at A 
and tracing out the corresponding z-plane behavior. 




1/3 scalar perturbation 

5/3 electromagnetic-gravitational perturbation 



(55) 




k < 1. 



The matching of asymptotic expansions in this case can be done following 
the recipe given in Section 2.5. As before, we begin at the point marked A 
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with the asymptotic behavior tp ~ (1,-1) a, and we have the boundary 
condition 

(-1)a = 0. (56) 
Going through a 2ir rotation in z gives (compare (32)) 

2ir : (m, n) -> (m + 4, n + 4) (57) 

which is in this case (1,-1)^ — > (5,3)^. Then we can go out along the 
lobe where in our approximation the behavior is purely oscillatory. After 
going around the lobe we have to return to the Bessel region a second time. 
However, the asymptotics in terms of z will be different this time: we chose 
the branch of z so that the first time we entered the Bessel region it was at 
z = 0, but now we have traveled an additional distance in z (from changing 
branches of the second logarithm in (51)), namely 

5 = -2mk 2 /{\ - k). (58) 

So we may write the solution as 

1p(z) = B + C + y / Uj(z - 5)Jj/ 2 (uj(z - 5)) + B-C-y/uj{z - S)J_j / 2 {oj{z - 8)). 

(59) 

Introducing a slight modification of our previous notation, we write 

(n 5 ) = ■ (n). (60) 

Then from (59) we obtain ip ~ (— lg, 1_s)b, giving the constraint 

(5,3)a = (-1*,1-«)b. (61) 

Making the second 2ir rotation we get (—1s,1-s)b —> (Ss,5-s)b, so the 
asymptotics near the point B are ip ~ (3<y, 5s)b- Now running over the big 
semicircle to get back to the point A, as before, the WKB approximation 
allows us to claim the coefficient of the exponentially big component e~ luJZ 
is unchanged. So over the full round trip this component is multiplied by 
(3<5)_b/(1)a- Between (56) and (61) we have three homogeneous linear con- 
straints on the four coefficients A±,B±, so we can determine this ratio; it 
turns out to be 



(1)a 



(l + 2cos^) -e (2 + 2cos7rj). (62) 



The contour encloses only the singularity at r = 1 (the outer horizon), so 
the monodromy comes only from the multivaluedness of ip and e~' luJZ around 
r = 1. As discussed in the previous section, this monodromy is determined 
by the surface gravity of the horizon at r = 1, namely 
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Also note that 2i5 = k 2 (3. Then matching the monodromies we obtain the 
final result, 

fP» = - (1 + 2 cos irj) - e k2 ^ (2 + 2 cos vrj) . (64) 

The first term looks exactly like the result (36) from Schwarzschild, and then 
there is a correction term which came from the phase shift 5. Actually, this 
term can also be written in terms of a Hawking temperature, but this time 
the Hawking temperature of the inner horizon: 

At any rate, substituting j = 1/3 or j = 5/3, from (55), gives the same 
answer in each case: 

e P» = -2- 3e fc2/ ^. (66) 

So this is our result for the asymptotic quasinormal frequencies of scalar or 
axial electromagnetic-gravitational perturbations of Reissner-Nordstr0m in 
d = 4. This result suggests a somewhat more complicated asymptotic be- 
havior for the quasinormal frequencies than we found in the case of Schwarz- 
schild; in particular, the asymptotic frequencies seem to be aperiodic unless 
k 2 is rational, and the real parts do not settle down to a simple value. 

We note that by taking k slightly imaginary it would be possible to get 
rid of the last term in the asymptotic limit, and in this case one would 
obtain asymptotically | Recj| = ^Hawking l°g(2). This may not be an unrea- 
sonable thing to do, since k is related to the charge and mass of the black 
hole, and the mass should be considered as slightly imaginary (the black 
hole eventually evaporates.) Indeed, if Q is real and fixed, then M is a 
decreasing function of k, so the expected small negative imaginary part of 
M corresponds to a small positive imaginary part of k, which is the cor- 
rect sign to suppress the last term. Following [1, 13] and identifying this 
^Hawking log(2) with a quantum of energy emitted by a black hole, we note 
that this corresponds precisely to the black hole reducing its outer horizon 
area by 4Glog(2) while keeping Q fixed. 

One might be concerned that our result for Reissner-Nordstr0m is incom- 
patible with the known results for Schwarzschild because it does not reduce 
to Tnawking log(3) as k — > 0. There is no inconsistency; rather, there is an 
order-of-limits issue. For extremely small k, Re 00 will be approximately 
^Hawking log(3) for a large range of Imw but then eventually switch over to 
the Reissner-Nordstr0m behavior determined by (66); the value of Imoj at 
which the crossover occurs goes to 00 as k — > 0. We can make a rough es- 
timate of the critical value k = ko(uj) above which the Reissner-Nordstr0m 
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description becomes more appropriate than the Schwarzschild description; 
namely, it is the point where the leading term fc 2 /r 6 is of the same order 
as the term 1/r 4 which would be leading for the Schwarzschild black hole. 
To compare these two terms the relevant r is the place where we patch the 
plane wave with the Bessel function solution, namely uz ~ 1. To relate z to 
r we can use either the relation z ~ r 2 from Schwarzschild or z ~ r 3 /k from 
Reissner-Nordstr0m. With either choice, requiring /c 2 /r 6 ~ 1/r , i.e. k ~ r, 
we obtain finally k ~ r ~ lj -1 / 2 i.e. 

jfco(w) ~ uj~ 1/2 . (67) 

If fc 3> fco(cj), the asymptotic Reissner-Nordstr0m formula (66) is appropri- 
ate; if k <C ko(uj) the Schwarzschild result is appropriate. 

After the preprint version of this paper appeared, the formula (66) as 
well as the order-of-limits issue discussed above were confirmed explicitly in 
the numerical analysis of [26] (see also [27] for some discussion.) 

4 Directions for further study 

This work suggests various possible directions for further study: 

• The machinery developed here may be applicable more generally, e.g. 
to charged black holes in d > 4, and perhaps even to general rotating 
black holes. Some work in this direction is in progress for the Kerr 
black hole in d = 4 [28]. In addition one could consider the case of 
black holes in de Sitter or anti-de Sitter space; quasinormal frequencies 
in these backgrounds have recently been studied [29, 30], and recently 
some exact results have been obtained in the limit where the Schwarz- 
schild horizon nearly coincides with the de Sitter horizon [31]. The 
asymptotic quasinormal frequencies in these backgrounds are related 
to properties of the holographic dual theory [32, 33]. The method will 
however require some modification, since the potential no longer van- 
ishes at infinity, and the expected answer is quite different (the real 
part is expected to diverge along with the imaginary part [29, 32].) 

• Our results so far only give the asymptotic real part of the frequency 
in the limit of infinite imaginary part. It may be possible to extend 
the technique to obtain a systematic expansion including corrections 
in l/y/\uj\. After the preprint version of this paper appeared, the j 
and I dependence of the first correction for arbitrary j appeared in 
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[27], and the exact correction for j = 2 including numerical coefficient 
was computed in [34]. 

• It would be desirable to have some numerical calculations of the asymp- 
totic quasinormal frequencies in the cases we considered above, to con- 
firm or reject our proposed analytic solution. In the case of Reissner- 
Nordstr0m, after the preprint version of this paper appeared, such 
calculations were indeed done [26] and confirmed our solution in detail 
(the agreement is also discussed in [27].) However, the case of higher- 
dimensional black holes remains numerically unstudied as far as we 
know. 

• In d > 4 Schwarzschild we could only give results for certain types of 
perturbation, because so far a good theory of more general perturba- 
tions in these dimensions is apparently lacking. It would be nice to 
understand how the leading term in the potential depends on the type 
of perturbation. 

• According to our calculation, the asymptotic real part of the frequency 
depends only on the coefficient of the leading term in the potential. 
This leading term becomes important only near the singularity at 
r = 0, but it still affects the transmission amplitude in the physical re- 
gion between the horizon and infinity. Although the importance of the 
singularity at r = was purely mathematical in our approach, there 
could be a more physical explanation of why our procedure works. Can 
we view it as a manifestation of the complementarity principle, reflect- 
ing some kind of duality between what happens inside and outside the 
horizon? 

• The occurrence of the Hawking temperature for the inner horizon of 
the Reissner-Nordstr0m black hole is surprising. How can the quasi- 
normal modes know about the behavior of the black hole in the causally 
disconnected interior? 

• The original motivation for this work was the LQG explanation [1] of 
the log (3) in d = 4 Schwarzschild in terms of a spin network with links 
of spin 1. If LQG is correct, perhaps the more complicated behavior for 
Reissner-Nordstr0m should have a similar explanation. The possible 
appearance of log(2) in the asymptotic frequencies could support some 
recent claims [35] that the gauge group of LQG should be SU(2) despite 
the log(3) for Schwarzschild. 

• Similarly, can spin foam models or other models explain the fact that 
the log(3) for Schwarzschild black holes seems to be independent of 
dimension? 
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We have answered some questions about the asymptotic quasinormal fre- 
quencies of black holes, but many more remain open. 
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